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Abstract In this paper, I present a general discussion of several astrophys- 
i-C ■ ical processes likely to play a role in the production of non-thermal emission 

in massive stars, with emphasis on massive binaries. Even though the discus- 
sion will start in the radio domain where the non-thermal emission was first 
detected, the census of physical processes involved in the non-thermal emis- 
sion from massive stars shows that many spectral domains are concerned, 
from the radio to the very high energies. 

First, the theoretical aspects of the non-thermal emission from early- type 
stars will be addressed. The main topics that will be discussed are respec- 
tively the physics of individual stellar winds and their interaction in binary 
systems, the acceleration of relativistic electrons, the magnetic field of mas- 
sive stars, and finally the non-thermal emission processes relevant to the case 
of massive stars. Second, this general qualitative discussion will be followed 
by a more quantitative one, devoted to the most probable scenario where non- 
thermal radio emitters are massive binaries. I will show how several stellar, 
wind and orbital parameters can be combined in order to make some semi- 
quantitative predictions on the high-energy counterpart to the non-thermal 
emission detected in the radio domain. 

These theoretical considerations will be followed by a census of results 
obtained so far, and related to this topic. These results concern the radio, 
the visible, the X-ray and the 7-ray domains. Prospects for the very high 
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energy 7-ray emission from massive stars will also be addressed. Two partic- 
ularly interesting examples - one O-type and one Wolf-Rayet binary - will 
be considered in details. Finally, strategies for future developments in this 
field will be discussed. 

Keywords Radiation mechanisms: non-thermal • Stars: early-type • Stars: 
binaries: general • Radio continuum: stars • X-rays: stars • Gamma rays: 
theory 



1 Introduction 

Stars with masses higher than about iO solar masses are often referred to as 
massive - or early-type - stars. This category includes O- and early B-type 
stars, along with Wolf-Rayet (WR) stars which are believed to be the evolved 
counterparts of O and B stars that have lost a substantial fraction of their 
mass through their strong stellar winds. 

Observational studies of massive stars reveal a wealth of crucial infor- 
mation concerning the physical processes at work in their expanding atmo- 
spheres. For instance, early-ty pe stars have been detecte d in the infrared and 
radio domains. Independentlv. lWright and Barlowl ( 19751 ) and lPanagia and Fellil 



( 1975( ) developed models to explain the production of the radio and infrared 



spectra of massive stars: the emission at these wavelengths consists of free- 
free radiation produced by electrons present in the plasma of their extended 
envelope. In this emission process, it is thus the stellar wind of the massive 
star - and not the star itself - that is the source of radiation. Moreover, 
the electrons involved in this process are distributed in energy according to a 
Maxwell-Boltzmann law: it is therefore a thermal emission process. The main 
characteristic of this thermal radiation is a continuum emission that can be 
described by a power law of the type: 

S v oc v a 

where S v is the flux density, v is the frequency, and a is the spectral in- 
dex of the free-free emission. For a homogeneous mass loss, at constant or 
accelerating velocity, a is about equal to 0.6, in agreement with many obser- 
vations. This thermal emission is intimately related to the mass loss rate of 
massive stars. A popular method to determine mass lo ss rates is moreover 
based on the measurement o f the thermal radio flux ( Wright and Barlowl 
119751 : iPanagia and Fellil Il975h . 

However, some of the massive stars observed in the radio domain present 
signifi c ant deviations from t his well-def i ned b ehaviour (see e.g. lWhite and Beckerl 
(|l983l ): lAbbott et al (|1984| ). or iRauwl (|2004 ) for a review). The main char- 



acteristics of such a radio emission disagreeing with the classical thermal 
emission discussed above are the following: 

- the spectral index (a) is not equal to 0.6, but is significantly lower than 
the thermal spectral index, and might even be negative. 

- the brightness temperature of the radio emission is much higher (~ 10 6 
- 10 7 K) than for thermal emission (~ 10 4 K). 
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- in most cases the radio emission is not steady, and the flux and/or the 
spectral index can present strong variations. 

This non-the rmal radiation was proposed to be synchrotron radiation by 
IWhitd (|!985h . Such a radiation mechanism requires two main ingredients. 



First, the existence of a magnetic field is needed. The direct detection of the 
magnetic field of massive stars is not an easy task, mostly because of the diffi- 
culty to investigate a Zeeman splitting in line profiles already broadened due 
to the fast rotation of the star. However, some recent observations led to the 
direct measurement of the magnetic field strength in the case of a few O- and 
B-type stars. Second, a population of relativistic electrons is required. The 
acceleration of particles in astrophysical sources has mainly been proposed 
to proceed through the Fermi mechanisn, in the presence of hydrodynamic 
shocks (see Section[3|). 

From this preliminary discussion, it becomes clear that one may discrim- 
inate between two main kinds of radio emission, produced by two distinct 
mechanisms whose physical requirements are very different: one is thermal, 
the other one is non-thermal. In this paper, a census of physical processes 
likely to lead to the production of non-thermal radiation will be made. The 
main aspects that I will address are respectively the physics of stellar winds 
and their interaction in massive binaries, the acceleration of particles, the 
magnetic field of massive stars, and finally the non-thermal emission pro- 
cesses likely to be at work in the context of early-type stars. 



2 Stellar winds and wind interactions 



Massive stars are known to produce strong stellar winds. In the case of O-type 
stars, these stellar winds can reach velocities up to 2000-3000 km s^ 1 with 
mass loss rates of about 10 -6 -10 -5 M Q yr . The mass loss from massive 
stars is driven by their strong radiation field, able to transfer momentum to 
the ions pres ent in their outer envelope . This line driving mechanism was first 
proposed bv lLucv and Solomon! ( 19701 ) and developed by I Cast or et all ( 19751 ) 
- the so-called CAK-theory - before be i ng further investigated and refined 
for instance bv I Abbott and Lucvl (fl~985h . IVink eTal (2000), P auldrach et all 
(|1986L 11994) andlKudritzkil |1996T) . A general discussion of this process can 
be found in De Beckerl (120011). 



In the basic CAK-theory, the stellar wind is con sidered to be stat i onary , 
i.e. the radial flow of matter is smooth. However, iLucv and Wiiitel (|l980l) 
pointed out the fact that radiatively driven stellar winds should be unsta- 
ble, leading to the appearance of perturbations in the mass outflow. These 
instabilities, generating intrinsic hydrodynamic shocks and clumps in stel- 
l ar wi n ds, were inve s tigate d in details for insta nce by lOwocki and RvbickH 
(|l985l) . IOwocki et all (|l988l ) and lFeldmeierl (|l995l ). The presence of the struc- 
tures resulting from outflow instabilities is believed to produce several obser- 
vational phenomena: 
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statistical fluctuati ons in emission li nes in the o ptical domain, in the case 
of O-type (see e.g . lEversberg et all (Il996l Il998l )) and WR stars (see e.g. 
iMoffat et all (|!994 : lBrown et all (|1993v T Such fluctuations have generally 



low amplitudes (of order a few percent at most) and individual 'clumps' 
affect rather narrow velocity intervals. 

the emission of thermal X-rays, likely produced by hydrodynamic shocks 
between plasma shells travelling at different velocities in the wind of iso- 
late d stars, and lead i ng to local temperature enhancements of a few 10 6 K 
(see iFeldmeier et al (|l997l )). The thermal nature of the X-rays produced 
within stellar winds is clearly revealed in their spectra by emission lines 
associated to highly ionized elements such as O, N, Si, Mg, Ne or Fe. 
Such emission line spectra ha ve been observed by the XMM-New ton and 
Chandra X-ray observatories (jSchulz et alll2000t Kahn et"all200lh . 



Beside these intrinsic radiatively driven instabilities, the winds of massive 
stars can be involved in strong collisions if the star is a member of a mas- 
sive binary system. In this scenario, the two companions produce their own 
stellar winds, that interact in a wind- wind collision zone somewhere between 
the two stars. This interaction zone is confined betw een two hydrodynamic 
shocks (jStevens et allll992t IPittard and Stevendll997Tj . 

A consequence of the colliding winds is the production of thermal X-rays, 
in addition to those produced intrinsically in individual stellar winds. The 
temperature of the plasma heated by this collision can reach a few 10 7 K, 
while that of individual stellar winds reaches at most a few 10 6 K. This 
additional emission component generally makes binary systems brighter in 
X-rays than the individual stars constituting them. Empirically, the ratio 
between the X-ray and bolometric luminosities of a single staiQ is about 
I0~ 7 , and any strong excess with respect to this ratio could be considered as 
an indication of an interaction between stellar winds in a binary system. 

In summary, as illustrated in Fig. [I] hydrodynamic shocks are common 
features in stellar wind physics, be they due to intrinsic instabilities or to 
wind-wind collisions in binary systems. This is a crucial point in the context of 
the non-thermal emission from massive stars, as it appears to be a mandatory 
ingredient required to accelerate particles. 



3 Particle acceleration mechanism 



When dealing with the hydrodynamics of stellar winds (individually, or in 
binary systems), an important concept worth introducing is that of strong 
shocks. Strong shocks are particular cases of hydrodynamic shocks where the 
pressure of the upstream gas (the gas that has not yet crossed the shock, 

1 A canonical relation between X-ray and bolometric lum inosities based on a larg e 
albeit heteroge nous samp l e of O -stars was proposed by !Bcrghocfcr ct al (1997). 
More recently, ISana etall lj2006T ) established a similar relation on the basis of a 
more reduced sample of stars that are members of the open cluster NGC 6231. The 
homogeneity of this latter sam ple led to a disp e rsion of the canonical relation of 
only about 40%, whilst that of Bcrghoefcr ct al (1997]) is about a factor 2. 
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Stellar wind 

Stellar photosphere Shock 




Fig. 1 Schematic view of the basics of the physics of stellar winds. The wind of a 
massive star is first produced through radiative driving. The stellar wind gives rise 
to instabilities that take the form of intrinsic hydrodynamic shocks. In the case of 
binary systems, the interaction of the stellar winds between the two companions 
can produce additional shocks in a collision zone. 



i.e. pre-shock gas) is negligible in comparison to that of the downstream gas 
(the gas that has already travelled through the shock, i.e. post-shock gas). 
Consequently, if the flow is considered to consist of a perfect gas, we obtain 
the following relation for the compression ratio fx)'- 

Pd v u 7 + 1 

X = — = — = 7 (1) 

Pu v d 7-1 

where pd (resp. Vd) and p u (resp. v u ) are the downstream and upstream gas 
densities (resp. velocities). For a monoatomic gas, the adiabatic coefficient 
(7) is equal to 5/3, and we derive one of the main properties of strong shocks, 
i.e. x = 4. 

A fi rst ap proach to explain the acceleration of particles was proposed by 



A nrst approacn to explain tnc acceleration 01 particles was proposed Dy 
iFermil (|1949I ). In his scenario, particles reflected between randomly moving 



"magnetic mirrors" (i.e. irregularities or turbulences in the magnetic field) 
are accelerated, resulting in a relative increase of their energy (Z\E/E) of 
the order of (V/c) 2 , where V is the velocity of the shock in the frame of the 
ambient plasma and c is the speed of light. 

A variant of this original idea, but in the pres ence of strong hydrody- 
namic shoc ks, was then developed for instance by iBlandford and Ostrikerl 
(|!978h and [Bel (jl978al lbh. In this context, the crucial point is that, from 



the particles' point of view, only head-on collisions occur and the accelerated 
particles do not undergo the decelerating effect of any following collision, 
contrary to what happens in the original scenario proposed by Fermi. In the 
literature, this acceleration process is also referred to as the Diffusive Shock 
Acceleration - DSA - mechanism. 

This is illustrated in Fig.fSJ where four situations are presented individ- 
ually. In situation (a), the upstream gas is considered to be stationary, and 
the shock crosses the medium at a supersonic velocity U. If the shock is con- 
sidered to be stationary (b), the upstream gas crosses it with a bulk velocity 
v u = U. However, as we are dealing with a strong shock (see Eq.Q]), the 
downstream bulk velocity Vd is equal to v u /4 = -j. In situation (c), i.e. 
from the point of view of the upstream gas considered to be stationary, the 
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downstream gas moves with a velocity v^l = U — ¥ — | {/. And in the 
reverse situation (d), i. e. when the dow nstrem gas is considered to be sta- 
tionary, v u ll = | U (see lLongairl ( 19921 ) for details). The main conclusion of 
this discussion is that the particles crossing the shock encounter gas moving 
with the same velocity, whatever the direction of the crossing of the shock. 
Consequently, the particle undergoes exactly the same process of receiving 
a small amount of energy A E when crossing the shock from downstream to 
upstream as it did in travelling from upstream to downstream. There are no 
crossings in which the particles lose energy : as a result, the process goes lin- 
early with V/c, and is therefore much more efficient than the original Fermi 
mechanism. 



a) 





U 


d 


u 



b) 



u 




d 


u 



c) d) 




Fig. 2 Illustration of the dynamics of high-energy particles in the presence of a 
strong shock where particles are accelerated through the first order Fermi mecha- 
nism. The upstream and downstream parts of the flow are respectively referred to by 
the letters u and d. The curves in panels (c) and (d) illustrate the scattering due to 
magneto-hydrodynamic phenomena, and responsible for the isotropic distribution 
of particles in the frame of reference where the gas is at rest. Th e four si t uation s 
are individually described in the text. This figure is inspired from Longair (1992). 



Another crucial aspect of this mechanism is that it is an iterative pro- 
cess. After travelling through the shock front, the high-energy particles can 
interact with magneto-hydrodynamic waves (Alfven waves) and magnetic 
turbulence. The consequence is an isotropic scattering of the particles which 
are allowed to cross the shock once again in the other direction. This results 
in the multiplication of the crossings, and therefore of the increases in en- 
ergy, which allows the particles to reach very high energies up to very large 
Lorentz factors. As a result of this iterative process, the energy spectrum of 
the relativistic particles is a power law: 

N(E) oc E- n (2) 

where n is the index of the power law. If we are dealing with electrons, 
this index will be referred to as the electron index of the relativistic electron 
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population. Following the study of lBelll ( 1978al ). the compression ratio of the 



shock and the index of the relativistic electron population can be related: 

X + 2 

n = _ (3) 

This general mechanism is believed to be responsible for the production 
of cosmic rays, probably at least partly accelerated in supernova remnants. 
In the context of stellar winds of massive stars, the first order Fermi mech- 
anism in the presence of hydrodynamic shocks (intrinsic or in wind-wind 
interaction zones) is a good candidate to explain the existence of the pop- 
ulation of relativistic electrons needed to produce the observed radio syn- 
chrotron emission. Figured summarizes the iterative process of acceleration 
of relativistic electrons in the presence of (i) hydrodynamic shocks and (ii) 
magneto-hydrodynamic phenomena likely to scatter isotropically the high- 
energy particles. This last point suggests that the role of the magnetic field 
of massive stars is not negligible in this acceleration process. 

Shock 

— Relativistic electrons 



Scattering by 
Multiple crossings . , „ 
Altven waves 



Fig. 3 Schematic representation of the Fermi acceleration mechanism in the pres- 
ence of strong shocks. The iterative nature of this process through the interaction 
of magneto-hydrodynamic phenomena is illustrated. 



Considering the role played by hydrodynamic shocks in the acceleration 
of particles, two scenarios can thus be envisaged. First, in the single star 
scenario, the relativistic electrons are accelerated by intrinsic shocks within 
stellar winds of isolated stars. Second, in the binary scenario, colliding winds 
are required to explain the non-thermal emission detected in the radio do- 
main. To evaluate which scenario is best adapted to explain the non-thermal 
emission from massive stars, the multiplicity of non-thermal radio emitters 
must be investigated. At this stage, we know that most of the non-thermal 
radio emitters - whatever their spectral type - are binaries or higher multi- 
plicity systems (see Tables [T] and [2]). The binary scenario appears therefore to 
be the most adapted to explain the non-thermal radio emission from massive 
stars. In this scenario, the non-thermal emission is produced in the wind- 
wind interaction zone between the two stars of a bi nary system, whe re the 
relativistic electrons are supposed to be accelerated (Dc Becker 2005). Re- 



cently, models have been developed applying the DSA mechanism to the 
partic ular physical conditions of hydrodynamic shocks of colliding-wind bi- 
naries ( Pittard and Doug hertvlbOOdrReirner et a]|l2006[ ). In the remainder of 



this paper, I will thus consider that we are dealing with massive binaries. 
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4 Magnetic fields of massive stars 

The magnetic field of stars is certainly one of the most poorly understood 
aspects of stellar physics. Observations of the Sun brought a huge amount 
of relevant information on stellar magnetism, but our general understanding 
of stellar magnetic fields is still very poor. However, as our knowledge of 
stellar astrophysics increases, the importance of the role played by magnetic 
phenomena becomes more and more obvious. In this context, it is worth 
considering in more detail the generation of magnetic fields within stellar 
interiors. 

Section 14.11 provides a short discussion of some fundamental hydromagnetic 
concepts. Sections 14.21 and 14.31 are devoted to hydrodynamic processes likely 
to be responsible for th e presence of a s tella r magnetic field. T hese discussions 
are mostly inspired by iTassoul ( 19761 ) and ICampbelll ( 19971 ). An alternative 



approach to produce a magnetic field is briefly described in Sect. 14.41 Sec- 
tion l4.5l summarizes the various processes possibly at work in the particular 
case of massive stars. Finally, Sect. 14.61 gives some words on observational 
methods used recently to estimate the strength of surface magnetic fields of 
a few massive stars. 



4.1 Hydromagnetic considerations 



Let us first consider some general equations of electromagnetism. For an 
electrically conducting fluid such as a plasma, Maxwell's equation relating 
the curl of the magnetic field (B) and the current density (J) is the following 



V x B — fi J 



(4) 



where /i is the permeability of free space. An additional equation useful for 
our purpose is Faraday's law 



V x E = 



dB 

~dt 



(5) 



If we include the induction electric field in Ohm's law, we obtain 

J = a(E + v x B) (6) 

where er, E and v stand respectively for the electric conductivity of the fluid, 
the intensity of the electric field, and the velocity of the flow. 
By taking the curl of Eq.[Hl we obtain 



V x J 



VxE + Vx(vxB) 



Using Eq.|3]in the left hand side and Eq.[5]for the first term of the right hand 
side yields 

V x B ^ 

V x (v x B) 



V x 
1 



Mo 



dB 

~dt 



dB 



-V x (V x B) = - — + V x (v x B) 
jU er ut 
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where we considered the permeability of free space as being constant. If we 
define the magnetic diffusivity as A = and use the vector identity 

V x (V x £?) = V • V • .B — V 2 £ where the first term of the second 
hand side is equal to zero because B is divergenceless, we finally obtain the 
induction equation: 

a d 

— = V x (v x B) + XV 2 B (7) 
at 

This equation expresses the temporal variation rate of the magnetic field 
in a plasma of finite conductivity through two terms: (i) the first one, also 
called the motion term, corresponding to a gain in intensity, and (ii) the 
second one stan ding for the dec ay of the magnetic field due to the resistivity 
of the medium (|Campbelllll997f ). 



4.2 Fossil magnetic field 

It has been proposed that the magnetism of some stars could be the relic 
of a magnetic field that was present in the g as cloud f rom w hich the star 
was formed by gravitational collapse (see e.g. iTassoull ( 19761) for a discus- 



sion). When trapped in a star without any mass motion, such a magnetic 
field would decay due to the finite conductivity of the stellar material. This 
can be illustrated by the diffusion equation, which expresses the decay of the 
magnetic field without any source of magnetic energy. 

If one neglects induction in Ohm's law, the first term on the right-hand 
side of Eq. [7] vanishes, and one obtains an estimate of the diffusion time-scale 
(jd) through the following relation 

I 2 

r d = j (8) 

where I is the typical length scale of the plasma. Consequently, if the plasma 
is turbulent, I will decrease substantially, leading to a rapid decay of the 
magnetic field. 

4.2.1 The Sun 



In the case of the Sun, ICowlind (|!945ft showed that the time-scale of electro- 
magnetic decay is about 10 10 years, i.e. the same order of magnitude as its 
evolution time-scale on the main-sequence. This suggests that any putative 
primeval magnetic field could possibly be responsible for its current magnetic 
activity. Consequently, the fossil magnetism has been considered as a viable 
candidate to explain the magnetism of the Sun, and other stars as well. 
However, it is not yet clear whether or not a large-scale primeval magnetic 
field could survive the pre-main-sequence phases. Moreover, in the particular 
case of the Sun, it is commonly assumed that the magnetic field is con- 
tinuously generated, and that all traces of primeval magnetism have long 
since vanished. The reasons are the following: (i) a turbulent convection in 
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the surface layers would cause the decay of the fossil field on a very short 
time-scale (typically a decade), and (ii) the polar solar magnetic field shows 
periodic reversals in contradiction with the existence of a stable primeval 
magnetising. For these reasons, an alternative theory is needed. This theory, 
i.e. the hydromagnetic dynamo, is briefly described in Section[4T3] 



4-2,2 Chemically peculiar (CP) stars 

Chemically peculiar (CP) stars are objects displaying atypical abundances, 
along with abundance inhomogeneities over their surface that may be due 
to the diffusion of elements under the influence of the stellar magnetic field. 
The fact that only about 10 % of the near-main sequence stars in the same 
stellar mass range as CP stars display strong magnetic fields ( the so-calle d 
'10% problem') could possibly be explained by a fossil origin (|Mossll200.^ . 
The amount of magnetic flux trapped during star formation is tentatively 
expected to vary substantially, with only the high flux tail of the distribution 
resulting in stars with a detected magnetic field on the main sequence. Mag- 
netic fields of a few 10 2 to several 10 4 G have been measured for chemically 
peculiar Ap and Bp stars ([Michaud et allll98ll : iHubrig et alll2004f l. 



4-2.3 White dwarfs and neutron stars 

The scenario of the fossil origin of the stellar magnetic field has recently 
been reconsidered in the context of the study of magnetic fields of compact 
degenerate stars such a s white dwarfs (W D) and neutron stars (NS). More 
detail s can be found in lTout et all ( 20041 ) and iFerrario and Wickramasinghel 



(120051) . 



In this context, a fossil magnetic field is believed to be trapped in stars 
during their formation process and the magnetic flux should be conserved. 
For stellar masses lower than about 2M©, the star becomes fully convec- 
tive during the Hayashi sequence and the fossil magnetic field has no chance 
to survive. This is e.g. the case of the Sun where the magnetism must be 
explained in another way (see the dynamo mechanism in Section l4~5]) . For 
masses larger than about 2M Q , the star does not experience a totally con- 
vective episode during the Hayashi phase and some fossil magnetic field may 
survive. The conservation of the magnetic flux as the stellar radius shrinks 
from that typical for main-sequence OB stars to that of a NS provides an 
apparently consistent relation between magnetic field strengths during the 
main-sequence (i.e. a few G to a few kG) and those observed for NSs (~ 10 11 - 
10 15 G). These numbers seem to lend some support to the conservation of the 
magnetic flux initially of fossil origin to explain the presence of the magnetic 
field in massive stars. In this scenario, inhomogeneities in the very weak mag- 
netic field present in extended regions of the Galaxy are likely to produce a 
large dispersion of stellar magnetic field strengths. 

2 We note that such a periodic reversal of the stellar magnetic field has also been 
obser ved in the case of another star, i.e. the K0 dwarf LQ Hydrae ( Kitchatino v et all 
[20001 1. 
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The conservation of the magnetic flux provides a possible interpretation 
to explain the presence of a magnetic field in WDs as well. In this ca se, 
magnetic fi eld strengths of abo ut a few kG ( Aznar Cuadrado et alll200"il ) to 
10 6 -10 9 G (|Schmidt et a]| l2003) have been observed. Such a dispersion in the 
magnetic field strength is also proposed to derive from heterogeneities of the 
local Galactic magnetic field. 



4-2.4 Massive non- degenerate stars 

As indicated above, the conservation of the primordial magnetic flux has 
been proposed to be responsib l e for the magnetism of massive stars (see e.g. 
iFerrario and Wickramasing hc (2005)). However, we can subdivide the inner 
structure of massive stars into two main zones: (i) the convective core where 
the fossil magnetic field is unlikely to survive, and (ii) the radiative envelope 
which results partly from accretion, and which is likely to harbour some 
residual fossil magnetic field. Moreover, the part of the massive star that 
shrinks to produce the NS is the core, and not the envelope. So, how can 
we explain that surface (main-sequence) magnetic fields fit so well core (NS) 
values after the supernova explosion? Is it fortuitous, or is there any strong 
physical relation between the convective core and the radiative envelope from 
the magnetic point of view? At this stage, this question is not answered, 
and this thus suggests that more than a simple magnetic flux conservation is 
needed to understand the processes responsible for the observed magnetism of 
massive stars, be they degenerate or not. Alternative scenarios are discussed 
in SectionsEO] and l4~4l 



4.3 Dynamo theory 

The classical principle of the dynamo action relies essentially on the fact that 
the motion of a conducting fluid (i.e. a plasma) is able to produce electro- 
motive forces capable of maintaining an initial magnetic field. The energy 
is supplied by the forces driving the fluid flow before being converted into 
magnetic energy. Consequently, as a first approximation, we will consider 
that the initial energy supply is huge, and mostly unaffected by this energy 
conversion process. 

In order to obtain the dynamo equation, we will use the same formalism as for 
the induction equation (Eq.[7]) . This dynamo equation will then be discussed 
in the framework of stellar magnetic fields. 



4-3.1 The dynamo equation 

For many astrophysical bodies, is shorter than their age, thus they need 
internal motions of conducting material to explain the observed magnetic 
fields that would otherwise have undergone ohmic decay down to negligible 
levels. This issue is investigated in the present section aiming at establish- 
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ing an equation similar to Eq.0 but specific to the case of the dynamo effect^. 

Let's consider an axisymmetric magnetic field that can be split into 
poloidal and toroidal components. According to 'Cowling's theorem', a steady 
axisy mmetric magnetic field cannot be maintained by dynamo action ( Cowling 
Il945l ). The problem is that in astrophysical bodies, the sym metry is dom i- 



nantly axial. A solution to solve this problem was proposed bv lParke"rl (1955), 



who noted that convection can cause rising cells of fluid to rotate and carry 
lines of force of the toroidal field into loops. Consequently, a large number 
of such loops could regenerate the poloidal component of the magnetic field. 
Such a turbulence is expected to arise from the Coriolis force, and is most 
commonly called the a- effect. 

The velocity and magnetic fields can be expressed as a sum of large-scale 
(mean field) and small-scale (fluctuating field) contributions. The a-effect 
is a small-scale turbulence acting on the mean (large-scale) magnetic field 
through a mean electromotive force (e) given by the following relation 

e = aB - /3V x B (9) 

where a is a coefficient that vanishes for homogeneous isotropic turbulence, 
and (3 can be regarded as an additional contribution to the fluid's magnetic 
diffusivity due to the turbulence (also called eddy diffusivity). 

Taking into account this additional electric field in the induction equation, 
we obtain the following equation for the mean field, i.e. the dynamo equation 

— = V x (v x B) + V x (aB) + (X + p)V 2 B 
at 

This equation allows axisymmetric solutions to exist for the mean magnetic 
field. As the magnetic diffusivity is essentially dominated by the eddy diffu- 
sivity, it is often neglected in the formulation of the dynamo effect. 
In this last relation, the magnetic and eddy diffusivities have been considered 
to be constant throughout the plasma. However, a more general formulation 
of the dynamo equation would be 



dB 

~dt 



= V x (v x B) + V x {aB) - V x (A + (3) V x B 



(10) 



With appropriate boundary conditions, this dynamo problem turns out to 
be an eigenvalue problem, aiming at finding eigenvalues (£) corresponding to 
non-decaying solutions of the type Z?(x,t) = i?(x) exp(£t), where x stands 
for th e spatial position vector. The study of ICharbonneau and MacGregorl 
(|200lh shows that the solutions for massive main-sequence stars (convective 
core + radiative envelope) point to magnetic fields that are concentrated in 
the immediate vicinity of the core-envelope interface. This is interpreted as a 
consequence of the fact that the magnetic diffusivity is expected to decrease 



3 This discussion is mostly inspired bv lTassouil (|1976l ) and ICampbelll (|1997D . 
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rapidly as one moves from the turbulent core to the stably stratified radiative 
envelope. 

The expression of the velocity can contain a term accounting for the 
differential rotation. In that case, the expressions for the poloidal and toroidal 
components of the magnetic field show that a magnetic field can be produced 
through two main mechanisms, respectively called the a- and i?-effect: 

a-effect. As stated above, it is the production of a magnetic field from 
turbulence arising from the Coriolis force. This process is able to produce 
both toroidal and poloidal fields. 

J?-effect. It consists in the production of a toriodal field through the 
shearing of the poloidal component of the field due to the non-uniform 
rotation of the star. 

The general case where both effects act simultaneously is often called the 
a 2 i?-dynamo. Particular dynamo cases can also be considered, such as the 
a 2 -dynamo where the J7-effect can be neglected, or the a^?-dynamo where 
the impact of the a-effect can be neglected with respect to that of the Q- 
effect in the production of a toroidal field from the poloidal field. 

Stars with very low mass (M < 0.4 M Q ) are essentially fully convective. 
The dynamo process can occur within the whole stellar interior. In the case 
of stars of solar and later types (M < IMq), the hydromagnetic dynamo is 
believed to occur somewhere within the outer convective zone. For instance, 
a magnetic field of a few Gauss is observed over most of the surfac e of the 
Sun, w ith strong localized fields of up to 4000 G in sunspots (see e.g. iTassoull 
(1976)). '. n this case, cyclic variations of the chromospheric activity suggest 
that the Sun-like dynamo produces a time-dependent field responsible for a 
cyclic magnetic activity. 

Fo r massive stars (M > 8-10 M Q ), the study of lCharbonneau and MacGregorl 
(2001) shows that such a dynamo action is possible in the core. However, the 
question is how a magnetic field generated by a core dynamo could travel up 
to the surface of the star where it can be detected directly (e.g. using spec- 
tropolarimetric methods, see Sect. 14.6]) or indirectly (e.g. synchrotron radio 
emission, putative magnetic stellar wind confinement). Two scenarios have 
currently been proposed to address this issue, and they are briefly discussed 
below. 



4-3.2 How to bring the magnetic field to the surface? 

The transport of the magnetic field of massive stars from the convective 
core to the surface constitutes a challenge for the stellar magnetism theory. 
Two scenarios have b een proposed: (i) the meridional ci rculation through 
Eddington-Sweet cells ( Charbonneau and MacGregorl200ll). and fii) the buoy - 



ant r ise of magnetic tubes along the rotation axis (jMacGregor and Cassinellil 
l2003h . 
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Meridional circulation - Eddington- Sweet cells. The rotation of early-type 
stars can lead to departures from sphericity, respons ible for pole-equator tem- 
peratu re differences in the radiative interior (see e.g. lCharbonneau and MacGregorl 
( 200ll )h These temperature differences cannot be balanced under the assump- 
tion of simultaneous radiative and hydrostatic equilibria. The reaction of the 
star is to generate a large-scale circulation in the meridional plane, known 
as the Eddington-Sweet circulation. As some of these flow streamlines come 
close to the surface of the star, and also close to the boundary between the 
convective core and the radiative envelope, Eddington-Sweet cells may be 
able to bring the magneti c field from the core up to the surface . 
According to the study of ICharbonneau and MacGregorl (l200lh . in the case 
of models with a strong magnetic diffusivity contrast between the core and 
the envelope (that are probably the most realistic ones), the dynamo effect is 
inhibited by the meridional circulation if it reaches a level suitable for a signif- 
icant magnetic flux transport. Moreover, such meridional circulation regimes 
are not likely to be attained in the interior of early-type main-sequence stars. 
Consequently, an alternative scenario is needed to carry magnetic fields from 
the stellar core to the surface. 



Buoyant rise o f magnetic t ubes. In this model of advection of magnetic flux, 
iMacGregor and Cassin elli (2003|) consider the dynamics of a circular ring, 
symmetric about the stellar rotation axis. This ring is an isolated concen- 
tration of magnetic field, which is allowed to begin an outward motion from 
the convective core thanks to the buoyancy. In this model, the ring, or mag- 
netic tube, is in mechanical equilibrium with the surrounding medium. We 
note however that the radiation pressure is neglected in such a simplified ap- 
proach, therefore restricting the discussion to stars with masses up to about 
10 Mq. The rotation is considered to be homogeneous, and no meridional 
circulation is taken into account. 

The time-scale of the rise of the magnetic tubes is shorter than the typi- 
cal evolution time-scale of massive stars on the main-sequence. However, the 
density deficit of such a ring is expected to decrease as it travels through the 
radiative envelope, causing the upward buoyant acceleration experienced by 
the ring to decrease as well. The ring indeed spends considerably more time 
(about a factor 100) travelling through a thin shell between 0.90 and 0.96 
R*, than in ascending from the core to a radius of 0.90 Ft,*. The inclusion of 
meridional circulation is expected to significantly decrease the time needed 
to transport the magnetic flux through the layers just below the surface. The 
joint action of the buoyancy transport and of the meridional circulation could 
therefore possibly lead to short rise times, able to transport a significant mag- 
netic flux to the surface during the very first stages of the main-sequence life 
of massive stars. 



4.4 An alternative dynamo scenario 



A completely different possibility worth considering is the production of a 
magnetic field through a dynamo mechanism just below the surface, without 
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the need of bringing it from the core. This approach is de scribed for in- 
stance bv lSpruitlj2002| ) or bv iMacDonald and MulTanl (|2004l ). showing that 



a dynamo can be driven by the differential rotation in a stable stratified zone. 



According to the study of lSp ruit (2002|), an instability of an initial toroidal 



field (even small) c reates a vert ical field component. Such an instability, called 
Tayler instability (|Spruitl [l999). consists of non-axisymmetric horizontal dis- 
placements of fluid elements. This vertical component of the field is distorted 
by differential stellar rotation which converts it into additional toroidal field 
lines. These horizontal field lines become progressively closer and denser, 
therefore producing a much stronger horizontal field. 



T his mechanism is usually re f ered t o as the Tayler-Spruit dynamo ( Spruitl 



120021 : iMaeder and Mevnel 1200.1 l2004h . In the context of massive stars, this 
scenario could be envisaged within the radiative envelope, therefore produc- 
ing a magnetic field close to the surface. 



4.5 Summary 

Considering the scenarios discussed in this section, we see that the origin of 
the magnetic field in massive stars is not yet well established. Still, we can 
make a census of the physical processes likely to play a role in this context. 

A fossil magnetic field is unlikely to survive in the turbulent regime of 
the core of the massive star, but a classical dynamo action might produce a 
magnetic field in that part of the star. Such a classical dynamo cannot be 
at work in the radiative envelope. However, an alternative dynamo driven 
by the differential rotation of the star may produce a magnetic field close 
to the surface. Moreover, the possibility to bring to the surface some mag- 
netic field produced in the core thanks to the combined action of meridional 
circulation and buoyancy cannot be completely rejected. We also note that 
a contribution of a fossil magnetic field might be present in the radiative 
envelope, as it results partly from the accretion of interstellar material that 
has not experienced a fully turbulent episode. 



4.6 Direct determination of magnetic field strengths of massive stars 

In the presence of a magnetic field, the degeneracy of atomic sublevels is 
lifted, splitting a given spectral line into several components (Zeeman effect). 
This effect has also the particularity of polarizing the light emitted through 
radiative transitions between sublevels affected by this splitting. The Zeeman 
effect is mostly efficient in circular polarization. Among the Stokes param- 
eters characterizing the state of polarization of quasi- monochromatic light, 
this circular component corresponds to the Stokes V parameter and responds 
to the line-of-sight component of the magnetic field (longitudinal Zeeman 
effect). It indicates the intensity difference between right- and left-handed 
circular polarizations. The other Stokes parameters respectively measure the 
unpolarized (I) and linearly polarized (Q and U) light. 
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In the case of massive stars, the strong rotational broadening due to the 
large projected rotational velocities (V rot sin i w 100 - 450 kins" 1 ) does not 
allow to detect directly the line splitting (or broadening) due to the Zee- 
man effect. Nevertheless, the Zeeman polarization can be observed through 
spectropolarimetric methods. Recently, on the basis of this principle, the in- 
tensity of the magnetic fields of only a few massive stars have been estimated. 
Three examples are respectively (3Cep (B1IV), 9 1 OriC (04-6V), and £Cas 
(B2IV). With this method, for a simple dipolar magnetic fiel d model, Bp 
and 6 for /3Cep are respectively of 360 ± 60 G and 85° ± 10° (jDonati et all 
l200lh. The same para meters for ^OriC are of 1100 ± 100 G and 42° ± 
6° dDonati et~aill2002lh and 335tg5° G and 18° ± 4° for (Cas dNeiner et all 
l2003h . In the case of these stars, more sophisticated models including higher 
order held components (quadrupolar,...) donot improve significantly the fit. 
Finally, we note that as this technique relies on a particular Zeeman com- 
ponent (the longitudinal one), it is sensitive to the phase of the rotational 
cycle of the star. Consequently, observing exposure times must be signifi- 
cantly shorter than the rotational period. For this reason, and also because 
of the low amplitude of the Stokes signature, this method is only applicable 
to bright sta r jj. In the fut ure, observations with spectropolarimeters suchs as 
ESPaDOns (|Donatill2004j ) mounted on the Canada-France-Hawaii Telescope 
(CFHT) are expected to improve significantly our knowledge of the strength 
of the magnetic field of massive stars. 



4.7 The magnetic field in the wind interaction zone of massive binaries 

If we consider that there is no magnetic field production in the collision 
zone of massive binaries, the local magnetic field - in the interaction zone 
- should come from the stars themselves. We need therefore to estimate the 
dependence of the magnetic field strength as a function of the distance to 
the star where it is produced. Some equations are given in Section [6~2l to deal 
with this dependence. In addition, we should consider that two stars may 
contribute to the local magnetic field and that the magnetic field may undergo 
some intensi fication in the post-shock plasma. This latter point has been 
addressed bvlBell (|2004h . who considered the amplification of the background 



magnetic field by the accelerated particles. Only a very simplified approach 
to estimate the strength of the magnetic fied in the colliding wind region is 
followed in SectionROl 



5 Non-thermal emission processes 

A crucial point in the context of this paper is the distinction between ther- 
mal and non-thermal emission processes. In thermal processes, the energy 
distribution of the electrons involved in the photon production mechanism 

4 The method was therefore not applicable to the cases of the stars discussed 
in Sect .[8] though a direct determination of their surface magnetic field strength 
would be crucial. 
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(e.g. thermal free-free emission, radiative recombination and collisional exci- 
tation followed by radiative decay) can be described by a Maxwell-Boltzmann 
law. In the case of non-thermal processes, particles (electrons, protons...) are 
relativistic and are distributed according to a power law (see Eq.[5]). 

Among the non-thermal emission processes, we can distinguish particle- 
field interactions (particles interacting with radiation or a magnetic field) 
and particle-matter interactions. The most common processes arc briefly de- 
scribed below. 



5.1 Particle-field interactions 
5.1.1 Synchrotron radiation 

A relativistic electron travelling in a magnetic field will be subjected to a 
Lorentz force orthogonal to the magnetic field direction. Consequently, the 
electron will adopt a helical path, and will emit photons {synchrotron radia- 
tion) . 

The characteristic energy of the photons is given by 



3 h e B sin 9 

Ec = 1 

4 7T m e c 



E 



m e c 2 



(11) 



where h is Planck's constant, e the electron charge, B the magnetic field 
strength, 9 the pitch angle, m e the mass of the electron , c the speed of light, 
and E the energy of the electron ( Cheng and Romero! 120041 ). This relation 



shows clearly that very energetic particles and strong magnetic fields are 
needed to produce synchrotron radiation in the very high energy domain. 
However, with rather modest magnetic fields and "moderate" relativistic 
electrons, synchrotron radiation can easily be produced in the radio domain. 
In the context of massive stars, the non-thermal radio emission discussed 
in the intro duction of t his work is most probably produced through this 
mechanism (|Whitdll985l) . 

It should b e noted that th e so-called Razin-Tsytovich effect (or Razin 
effect, see e.g. IVan Lool (|2005l )) is likely to occur, leading to an apparent 
suppression of the synchrotron radiation. This effect is based on a particu- 
lar property of the synchrotron emission: it is characterized by a significant 
beaming, i.e. the radiation is emitted in preferential directions. This prop- 
erty of the synchrotron emission is important considering the fact that, in 
a plasma such as that constituting stellar winds, the beaming effect can 
be substantially reduced, leading therefore to a decrease of the synchrotron 
emitting power in a given direction. This suppression is a characteristic ef- 
fect of plasma on radiation by relativistic particles. It can be shown that the 
Razin-Tsytovich effect is a consequence of the suppression of radiation with 
a frequency below the plasma frequency, after transformation from the rest 
frame of the relativistic particle to the rest frame of the observer . A detailed 
discussion on this transformation can be shown in iMelrosd (| 19721) . This effect 



may be responsible for significant underestimates of the radio synchrotron 
power of astronomical sources. 
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5.1.2 Inverse Compton scattering 

Particles (electrons or protons) can also interact with photons. In the con- 
text of this study, only electrons will be considered. If the typical energy of 
the relativistic electrons is high enough, the latter can transfer some of their 
energy to soft photons. This process is known as inverse Compton (IC) scat- 
tering. In this way, moderately relativistic electrons are able to "promote" 
ultra-violet photons to the high-energy domain (X-rays and soft 7-rays). The 
relation between the average energy of the emerging photons {Eic) and that 
of the scattered soft photons (E so f t ) is 

Eic = 5 E soft 7 2 (12) 

where 7 is the Lorentz factor characterizing the relativistic electrons. This 
latter relation is only valid in the Thomson limit. In this regime, the square 
root of the product of the energies of the interacting ele ctron and photon is 
signifi cantly lower than the rest energy of the electron ( Cheng and Romero! 
12004 ). 

Most frequently, inverse Compton scattering is considered in the Thomson 
regime, and in that case the cross-section of the process is the Thomson cross- 
section (<tt, see Sect. 16 .5]) . However, when we are dealing with very high- 
energy particles, a quantum relativistic cross-secti on has to be u sed and the 
pr ocess is considered i n the Klein-Nishina regime ( Longairlll992f ). As shown 



by iReimer et al (2006), the Klein-Nishina cross-section can be significantly 



lower than the Thomson cross-section in the very high-energy domain. The 
main effect is a strong decrease in the efficiency of IC scattering to produce 
very high-energy photons. 

IC scattering is of particular interest in the framework of massive stars 
because of the strong supply of UV and visible photons provided by the 
photosphere. This process is indeed believed to play a significant role in the 
cooling of relativistic electrons in stellar wind environments. 



5.2 Particle-matter interactions 
5.2.1 Relativistic Bremsstrahlung 

When a relativistic electron is accelerated in the electric field of a nucleus or 
other charged particles, radiation is produced: relativistic Bremsstrahlung. In 
general, only the emission from the electron has to be considered. It comes 
from the fact that the intensity of Bremsstrahlung by a nucleus of mass M 
is ~ (M/m e )~ 2 times the effect generated by the electron. 

We note that such a mechanism requires rather high plasma densities 
to become significant. For this reason, relativistic Bremsstrahlung is not ex- 
pected to be responsible for non-thermal emiss ion in th e case of O-type stars. 
This scenario was nevertheless considered by iPollockl ( 19871 ) in the case of 
the denser winds of Wolf-Rayet stars. 
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5.2.2 ir° -decays from proton-proton emission 

Even though we are mostly dealing with electrons in the context of this pa- 
per, we should not neglect the fact that protons are very abundant in astro- 
physical environments such as stellar winds. Relativistic protons can interact 
with other 'ordinary' protons to produce neutral pions (ir ), that decay into 
gamma-rays with a proper lifetime of only 9 x 10~ 17 s. For a more detailed 
discussion of this process I refer to [Ch eng and Romero! ((2004). 

It has been propose d that the 7r°-decay m echanism could be at work 
in young open clusters ( Manchanda et all Il996h . In this context, it is sug- 



gested that the combined winds of massive stars may interact with the gas in 
the cluster and produce a system of bow shocks. These shocks could then be 
responsible for the acceleration of protons up to relativistic velocities. Conse- 
quently, the relativistic protons may interact to produce a diffuse gamma-ray 
emission in open clusters and OB associations. In addition, it is important 
to note that the 7r°-decay mechanism may certainly be responsible for the 
very high-energy emission discussed later in this review. 



6 The non-thermal high-energy counterpart of a massive binary 

I have discussed several aspects of the physics of massive stars, or differ- 
ent topics of astrophysics, likely to play a role in the production of non- 
thermal emission of massive stars. Figure|4] presents an overview of these 
physical processes. This schematic view is mostly inspired from that given by 
iDe Becker et"a1 (l2005bh . In this figure, the four boxes individually represent 



the topics independently discussed in the four previous sections of this chap- 
ter in their order of appearance. The intersections between boxes illustrate 
the noticeable interconnections between these different fields of astrophysics. 

As shown in Sect. El non-thermal emission processes can also be at work in 
the high-energy domain. In order to investigate the issue of the high-energy 
counterpart to the non-thermal emission of massive stars detected in the ra- 
dio domain, it is helpful to adopt a somewhat more quantitative approach of 
the physical processes believed to be responsible for the non-thermal emis- 
sion of radiation in the particular case of binary systems. We note that the 
description below should be considered as a simplified approach, used only to 
give an idea of the impact of several physical parameters on some quantities 
such as the IC luminosity. A more complete and quantitative approach would 
require sophisticated models including more physics than used below. 

The radio synchrotron luminosity and the inverse Compton luminosity 
can be related by a unique equation, provided both emis sion mechanisms 
arise from the same population of relativistic electrons (jWhite and Cher] 
I1995D . Such a relation derives from the fact that, for both processes, the 
total energy lost by a high-energy electron can be expressed by similar rela- 
tions depending on the energy density respectively of the magnetic field (Ub) 
and of the radiation field (U p h). This leads to the following equation for the 
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Fig. 4 Overview of the physical processes likely to contribute to the production 
of non-thermal emission from massive stars. The main aspects of the problem are 
subdivided in boxes directly related to the processes discussed with more details in 
this paper: (1) stellar wind physics (Sect.f2}; (2) particle acceleration mechanism 
(Sect.[3]); magnetic processes (Sect. [4}; non-thermal emission processes (Sect.[5]). 



ratio between the radio synchrotron luminosity (L sync h) and the high-energy 
inverse Compton luminosity (Lie)'- 

L synch Ub (13) 



Lie U p h 

We can use this luminosity ratio to obtain an expression of the IC luminosity 
as a function of the synchrotron radio luminosity, taking into account the 
fact that both stars of the binary system contribute to the photon energy 
density: 

L ic = -^^{— + —) (14) 

where Lb ^i (respectively Lbd^) and r\ (resp. r-i) stand for the bolometric 
luminosity and for the distance from the center of the star to the collision 
zone for the primary (resp. the secondary). B co i is the strength of the mag- 
netic field in the acceleration region of relativistic particles (which is also the 
emitting region of the synchrotron and IC ph otons, i.e. the colli s ion zo ne). 

This approach is very similar to that of lEichler and Usovl (jl993h and 
Bcna glia et all (2001b). However, these latter studies took into account only 



one term for the photon density. This approximation was valid because the 
authors had to deal with asymmetric binary systems where the winds collide 
close to the secondary. Consequently, the primary photon density could be 

general 
and 



neglected compare d to that of the secondar y. Eq.fMlis therefore more ge 
than that used by lEichler and Usovl (|l993h and Benaglia et aj (2001b) 



it should preferentially be used each time one has to deal with binary systems 
characterized by a momentum ratio (rj) that is not significantly smaller than 
1 (see Eq.ri9|) . It should be mentioned that the simplified approach described 
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below does not take into account the variation of the Us/Uph ratio accross the 
wind collision region, even though we know that the non-thermal emission 
region can be extended (see Sect. l8.ip . A more consistent approach would 
require a detailed knowledge of the spatial distribution of relativistic particles 
- along with that of the magnetic field - in the wind collision region. In the 
remaining of this section, I will discuss in detail three points that need to be 
addressed properly if one wants to apply Eq.[l4l the geometry of the system, 
the intensity of the local magnetic field, and the synchrotron radio luminosity. 



6.1 Geometry of the system 

6.1.1 Separation between the stars. 



In an eccentric binary system, the relative distance between the two stars in 
units of the semi-major axis of the orbit, is given by 



1-e 2 



1 + e cos v 



(15) 



where v is the true anomaly, i.e. the angle between the vector radius and 
the direction to periastron, and e is the eccentricity of the orbit. In order 
to obtain the relative separation as a function of the orbital phase (</)), a 
relation between the true anomaly and the orbital phase is needed. As a first 
step, let us consider the expression of the true anomaly (v) as a function of 
the eccentric anomaly (E). 



= 2 arctan 



tan 



E 



The eccentric anomaly is related to the mean anomaly {M = 2 n < 
= at periastron) through Kepler's equation: 



(16) 
■ when 



E - e sin E = M (17) 

The numerical resolution of Eq.[T7] allows to obtain the value of E for a 
given orbital phase. Then the successive application of Eqs. 1161 and 1151 pro- 
vides the relative separation between the two components of the binary sys- 
tem. 

The projected semi-major axes (ai sin i and a-i sin i) for the two compo- 
nents can be obtained from the spectroscopic orbital solution of the binary 
system. Provided an estimate of the inclination angle i is available, the multi- 
plication of the relative separation (r) by the sum a± + yields the absolute 
separation (D) between the two stars. 



6.1.2 Wind momentum equilibrium. 



The collision between the stellar winds of the two components of a binary 
system occurs at a position which is set by the relative strengths of the winds. 
If the winds have similar strengths, the collision zone is a plane perpendicular 
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to the line of centres, at equal distance from the two stars. On the contrary, 
if the winds are not equal, the collision zone is a curved layer of shocked gas 
partially folded around the star producing the weaker wind. 

If we assume that the spherical winds of the two components flow nearly 
radially out to the shocks, we can derive the distances from the stars where 
these winds meet. At the position of the collision zone, the wind momenta 
are balanced, yielding 

n = —L-rD,^ = -^tD (18) 

where r± , r-i , D and 77 are respectively the distances to the collision zone from 
the primary and the secondary, the separation between the stars (see above), 
and a dimensionless parameter expressed as 

77 = —. (19) 

Mi V x ,x 

where M and stand for the mass loss rate and the terminal velocity 
respectively. The parameter 77 thus ex presses the ratio of the momentum flux 
of the two winds ( Stevens et a 



6.2 Local magnetic field. 



The surface magnetic field strength (B s ) is probably one of the most critical 
and poorly known physical quantities for massive stars. Moreover, provided 
that B s is determined, the behaviour of the magnetic field as a function of 
the distance to the star in the ionized material of the stellar wind still needs 
to be established. 

Some relations providing the dependence of the strength and geometry 
of the magnetic field as a function of the radia l distance (r) to the sta r 
in the outflowing gas are given for instance by lUsov and Melrose! (|l992h . 
Three regimes can be envisaged: (1) dipolar, (2) radial, and (3) toroidal. The 
equations for these three regimes, along with their conditions of application, 
are given below: 



B = B, x < 



for R < r < r a (dipolar) 
^3 forr A < r < R^ (radial) 
< r (toroidal) 



(20) 



v^RL forR^- 

Vrot r A r V rot 



In these relations, V ro t and Vx, are respectively the stellar rotational 
velocity and the wind terminal velocity. R is the stellar radius, and ta is the 
Alfven radius defined as the radius where the kinetic and magnetic energy 
densities are balanced. The Alfven radius can be expressed as follows: 
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rA ~ iJ x \ e/i >> j (21) 



with the parameter £ expressed by 

B 2 i? 2 



2MK, 



(22) 



From these relations, we can see that the typic al Alfven radius for early -type 
stars is in the range from ~ i? to ~ 2 — 3 i? (Usov and Melrose 1992). The 



best approach consists in using the well-established or expected wind and 
stellar parameters to estimate which regime is suitable, before applying the 
corresponding relation to estimate the local magnetic field strength at a given 
distance r from the center of the star. 

When the strength of the magnetic field at the distance of the shock zone 
is obtained from Eq.J^Ol one still has to wonder whether it is affected by the 
hydrodynamic conditions close to the shock. In the case of a strong shock, the 
compression ratio (%) is expected to be equal to 4 (see Sect. [2]). The magnetic 
field strength dependence on the density could be of the order of p 2 ^ 3 ( Tassoull 



Il976l ). and can be used to determine some kind of concentration factor e. 
Consequently, the magnetic field strength could increase by up to a factor 
of 2.5 in the post-shock region of a strong shock. Moreover, the magnetic 
fields from both stars are expected to contribute to the local magnetic field 
strength at the collision zone. Let us therefore consider that the mean local 
magnetic field can be obtained using the following relation: 

B col = e^±^l (23) 

where B\ and B 2 are the local magnetic fields from the primary and the sec- 
ondary respectively, e xpressed in Gauss. This a ppro ach is thus sign i ficantly 



onciary respectively, e xpressed m uauss. inis a ppro acn is tnus sign ificantly 
different from that of lEichler and Usovl (| 19931) and iBenaglia etaJ (|2001bf ) 



where the local magnetic field was considered to come from the secondary 
star only, as those authors dealt with cases where the collision zone was very 
close to the secondary. I emphasize on the fact that the way we deal here 
with the magnetic field is only an assumption. A detailed investigation of 
the magneto-hydrodynamic conditions at the collision zone is far beyond the 
scope of the present discussion. 



6.3 Radio luminosity. 

In the literature, radio emission levels are mostly provided as flux densities 
expressed in mJy at a given frequency (S v ). This flux density can be con- 
verted into more conventional units (cgs) following the procedure described 
below. 

We first need to establish a relation yielding the flux in a given frequency 
interval. This can be done through an integration of the flux density between 
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two frequencies v\ and z^: 




S v dv 



(24) 



The lower boundary v\ is denned by the rather steep increase of the radio 
flux as a function of frequency likely to occur at frequencies of the order of 
a fraction of GHz (i.e. at a wavelength of a few tens of cm), before reaching 
a maximum at a frequency of the order of the GHz. This turnover may be 
due to the combined action of the opacity of the plasma that decreases with 
increasing frequency, and of the Razin-Tsytovich effect. Recent modelling 
and observations seem to fa yor free-free opacity as the dominant factor at 
least in the case of WR 140 Becker et a 1 2001E). For the upper limit on 
the frequency (z^), a value of the order of several tens of GHz ( i .e. a f raction 
of cm) should be considered (see e.g. Sect. 2.3.1 in lPittard et all (|2006f) ). This 
upper boundary is determined by the maximum energy reached by relativistic 
electrons that is severely affected by IC cooling (maximum Lorentz factors 
for electrons of the order of 10 4 -10 5 are expected depending of the colliding- 
wind binary system, see Sect. 16. 5[) . It should be noted that these characteristic 
frequencies depend strongly on the target that is considered. A few examples 
of synchro tron radio spectra of long perio d m assive binaries can be found for 
instance in lPittard and Dougherty! (|2006r ) and lDoughertv et all (|2003T ). As we 
know that the radio spectrum is a power law with a spectral index a, we can 
use the following relation: 



where v' is a frequency for which we have a measure of the flux density. 
Inserting this relation in Eq.[Ml we obtain 



As a final step, we consider that 

lmJy = 10" 29 Wm- 2 Hz _1 = 1CT 26 ergs _1 cm" 2 Hz _1 

and we multiply both sides of Eq.[5H] by 4 7r d 2 in order to obtain a relation 
giving the radio luminosity in cgs units integrated between two frequencies 
V\ and V2'. 





and we derive the integrated flux: 




(25) 
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v. 



■ la 



where S v >, v' and d are respectively expressed in mJy, in Hz and in cm. To 
calculate this quantity, the flux density at a given frequency (z/), the distance 
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to the star, and the spectral index are thus needed. 

In a binary system producing non-thermal radio emission, the total ra- 
dio spectrum has essentially three origins: (1) the thermal emission from the 
wind of the primary, i.e (2) the thermal emission from the wind of the 

secondary, i.e £//,2, and (3) the non-thermal emission from the collision zone 
between the two stars, i.e L sync h- This third part is likely to be embedded 
inside the winds and therefore undergoes a free-free absorption characterized 
by a given optical depth (r) . We note that a thermal ra dio emission compo - 
nent is also expected to be produced by colliding winds ( Pittard et all l2006). 



but we neglect it in our simplified approach. Consequently, the observed radio 
luminosity can be expressed as follows: 

L ° h adio = L ffA + L ff,2 + L synch e~ T (27) 

In order to estimate the non-thermal contribution to the radio spectrum, 
we need to estimate the free-free thermal contributions from the primary and 
the secondary respectively. These q uantities can be calcu lated on the basis 
of the relations given for example bv lLeitherer et all (11995), according t o the 



free-free radio emission t heory established by I Wright and Barlowl ( 19751 ) and 
iPanagia and Fellil (|l975l) : 



where g v is the Gaunt factor obtained through the following relation: 



g v = 9.77 1 + 0.13 log 



rpZ/2 

Zv 

The various quantities used in these two relations are 

M : the mass loss rate in M Q yr" 1 
Vqo : the terminal velocity in kms -1 
v : the frequency in Hz 
d : the distance to the star in kpc 

T e : the electron temperature, which is generally a fraction of the effective 
temperature, e.g. T e = 0.5 x T e ff 

H : the mean molecular weight, i.e. ^=f=^ - — -, where Ai and Mi are respec- 

tively the abundance and the molecular weight of the ith ionic species. 
For a plasma containing 77% hydrogen and 23% helium, fi — 1.69 

(V At Z 2 ) 1/2 

Z : the rms ionic charge, i.e. \ l f , where Zi is the ionic charge of 

the zth ionic species. For the same abundances as above, Z =1.3 

Ai Zi 

7 : the mean number of electrons per ion, i.e. ^ ' . For the same 
abundances as above, 7= 1.23 

Following this approach, the thermal flux density from both components 
of the binary system can be calculated. Using Eq. 1261 with a spectral index a 
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of 0.6, these quantities can be converted into radio luminosities expressed in 
ergs -1 . 

Applying Eq.[26] using S v i and a obtained from the observation, the ob- 
served radio luminosity can be calculated as well. The difference between 
these observed (mixed thermal and non-thermal) and theoretical (thermal) 
quantities provides a lower limit on the luminosity {L sync h y min) of the intrin- 
sic synchrotron emission produced at the collision zone: 

Lsynch,min — L synch & 

where r is an optical depth characteristic of the total absorption undergone 
by the synchrotron radiation produced globally in the emission region lo- 
cated at the collision zone. This optical depth results from the integration 
of the opacity along the line of sight over the entire emitting volume. The 
determination of such a quantity is not straightforward and requires an accu- 
rate knowledge of the geometrical and wind parameters of the system. This 
absorption factor (e~ r ) is likely to have a large impact on the synchrotron 
luminosity, mostly in the cases of close binary systems and systems including 
stars with dense stellar winds (e.g. WR stars). For instance, an optical depth 
t of about 2.3 would be responsible for an attenuation of the synchrotron 
luminosity by a factor 10. Moreover, we note that since we neglect here the 
Razin effect (see Sect. l5.lTT|) the intrinsic synchrotron luminosity will be un- 
derestimated. 



(29) 



6.4 Estimate of the lower limit on the inverse Compton luminosity 

Provided that (1) the bolometric luminosities of the stars of the binary sys- 
tem are known, (2) the distances between the stars and the collision zone 
have been calculated, and (3) the lower limit on the intrinsic synchrotron 
luminosity has been estimated, Eq.[l4] can be applied. It should be kept in 
mind that the intrinsic synchrotron luminosity is expected to vary with the 
orbital phase if the binary system is eccentric. This comes from the fact that 
the emitted synchrotron power depends on the strength of the local magnetic 
field and on the physical extension of the synchrotron emitting region, both of 
them varying with the orbital phase. When using the minimum synchrotron 
luminosity as defined by Eq.[29]at a given orbital phase, this approach leads 
to a lower limit on the high-energy IC luminosity at the same orbital phase. 



6.5 Estimate of the high-energy cutoff of the non-thermal spectrum 

The electrons are accelerated through the first-order Fermi mechanism, and 
a valuable quantity worth to be estimated is the maximum energy reached by 
the electrons under given circumstances. As non-thermal radiation processes 
provide the most effective cooling mechanisms of relativistic electrons, a high- 
energy cutoff is obtained when the rate of energy gain of the particles is 
balanced by the rate of energy loss by radiative processes like IC scattering 
and synchrotron emission. I note that relativistic electrons may also lose 
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energy through Coulomb collisio ns with thermal ion s, but this cooling process 
is not expected to be significant (jPittard et al 2006) : I thus decided to neglect 
it. 

On the one hand, we can use the relation given by lDoughertv et al ( 2003) 



for the rate of IC loss of energy of relativistic electro ns, e.g. electrons whose 
Loren tz factor 7 is 3> 1, and the relation provided bv lRvbicki and Lightmar] 
for the energy loss by synchrotron emission: 



dE 
~dt 

dE 
~dt 



ic 



<tt 7 
3tt 



<tt 7 



synch 



6 7T 



Lboi 
cB 2 



(30) 



(31) 



where <jt is the Thomson cross section whose value is 6.6524 x 10~ 25 cm~ 2 . 
Lboi and r are respectively the stellar bolometric luminosity and the distance 
from the star to the shock region. To take into account the radiation from 
the two stars of the binary system, the last factor of Eq.[30l should be written 



(- 



~ ) • 



On the other han d, the rate of energy gain through the Fermi acceleration 
mechanism given by iBaring et all ( 19991 ) can be expressed in cgs units as 
follows: 

%BVX (32) 



dE 
~dt 



A 



Fermi 



x(i + sx) 



where A is a constant ~ 1.6 x 10~ 20 . 

In this relation, \ is the compression ratio of the shock, B is the local 
magnetic field in Gauss, V s h is the speed of the shock in cms , g is the 
ratio between the upstream and downstream diffusion coefficients of the high- 
energy particles generally considered equal to 1/x, and 9 is the ratio of the 
mean free path of the particles over the gyroradius. This latter parameter 
is considered to be constant across the acceleration zone, to be independent 
of the particle and energy, and to characterize the turbulence. Its minimum 
value is one (for the case o f strong hydromagnetic turbulence), and it was 
considered to be about 3 bv lEichler and Usovl ( 19931 ) in the case of the wind 
collision zone of WR + O binary systems. A low 9 value means that the 
turbulence is high enough to reduce substantially the mean free path of the 
particles, therefore confining them in the acceleration region. Q is the mean 
electric charge. 

The particles are accelerated as far as the rate of energy gain is not 
overwhelmed by the rate of energy loss. The maximum Lorentz factor of the 
relativistic electrons is consequently provided by the condition: 



dE 
It 



ic 



dE 
~dt 



dE 

synch dt 



Fermi 



which translates to 



< 



3irA( X -l)QBV s \ 



^TX(l + .9X)(?(^f + ^f + 



(33) 
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The maximum Lorentz factor ("j m ax) expressed by the right hand side of 
Eq.[33] corresponds to an electron energy of 7 m e c 2 , where m e is the mass of 
the electron. Lorentz factors of the order of 10 4 -10 5 can be reached depend- 
ing on the long period binary system, with longer periods leading to higher 
Lorentz factors. For a -f max of the order or 10 4 , the maximum energy reached 
by the electrons is about 5 GeV. When the maximum Lorentz factor is esti- 
mated, one can determine the maximum frequency of the photons produced 
by inverse Compton scattering, through the following relation (see Eq. [T2")) : 

4 2 

V max = ^lmax V * ( 34 ) 

where is the typical frequency of the seed stellar photons (seerBlumcnt hal and Gouldl 
(|l970t) 1. This quantity can be estimated from Wien's law providing the fre- 
quency of the maximum of Planck's function for a black body temperature 
corresponding to the effective temperature of the star. For massive stars with 
a typical temperature of about 40000 K, the radiation field peaks at about 
4 x 10 15 Hz, leading to a maximum energy for the Comptonized photons of 
the order of a few GeV. 

In Eq. 1331 both inverse Compton scattering and synchrotron emission are 
considered. However, it should be noted that in most cases only IC scattering 
will significantly contribute to the energy loss of the relativistic particles, 
and the third term between brackets in the denominator of Eq.[33] can be 
neglected. The reason is that the UV and visible radiation fields from the 
stars are strong, and that the synchrotron mechanism would require too high 
a local magnetic field strength to become significant. 



7 The multiplicity of non-thermal radio emitters 

As mentioned in Sect.GU the binary scenario appears to be the most adapted 
to explain the non-thermal radio emission from massive stars. This statement 
comes mainly from the fact that the binary fraction among these particular 
sources of radiation is significantly large. 

First, let us consider W R stars. A first census of the non -thermal emitting 
WR stars was presented bv lDoughertv and Williams! ( 20001 ). who pointed out 



the fact that a large fraction of them were binary systems. At the time of the 
writing of this paper, 17 WR stars have shown a non-thermal contribution in 
their radio spectrum (see Table[T]). In most cases, the non-thermal nature of 
the radio emission is revealed by the spectral index that deviates significantly 
from 0.6. In some cases, the non-thermal emission is strongly suggested by 
significant variations in the radio flux. A striking feature of the list provided 
in Table[T]is the large fraction of double or multiple systems among the non- 
thermally emitting WR stars. Out of 17 non-thermal radio emitters, 14 are 
at least suspected binaries, with most of them (12) being confirmed binaries. 

In the case of O-type stars, the situation is rather similar. Up to now, 
16 non-thermal radio emitters of spectral type O are known (see Table[2]). In 
this category, 11 are confirmed binaries and two more are suspected binaries. 
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However, as the situation is less clear than for WR stars, the references for 
their mutiplicity is also given in TableO 

In order to evaluate the significance of the binary scenario for non-thermal 
radio emitters, it is useful to give a few words on the binary fraction of mas- 
sive stars in general. As most massive stars are members of open clusters, 
an estimate of their binary fraction can be obtained through multiplicity 
st udies of massive stars in open clusters. Some results have been published 
by iGarcia and Mermilliodl ( 200lh . resulting in intra-cluster binary fractions 
ranging between 20 and 80% depending on the open cluster. However, in- 
depth studies of the clusters with the largest cl aimed binary fractions showed 
that values of 80% are heavily overestimated. ISanal (l2007f) derived indeed a 
binary fraction of the order of 60 % in NGC 623 1 , and lDe Becker et all (|2006ah 
showed that the binary fraction among massive stars in IC 1805 should range 
between 20% and 60%. Putting these numbers with results related to other 
clusters, the binary fraction among massive stars in general should not be 
higher than 50 % - and this number should be considered as a very conserva- 
tive upper limit. The lower limit on the binary fraction among non-thermal 
radio emitting massive stars is therefore significantly higher than the upper 
limit on the binary fraction of massive stars in general. Studies are however 
currently being carried out in order address in detail the issue of the mul- 
tiplicity of massive stars in general, and of non-thermal radio emitters in 
particular. 



Table 1 List of the non-thermal radio emitting WR st ars as of November 2005 
The in for mation collected i n this tab le is taken from (a) Dougherty and William! 
pJOOh , fbllChapman et all (Tl999l ). fc) IBenaglia"et~a] (|2005l ), (d) ICappa et all | 200~ 
and (e) IMonniereTa l (2002). 



Star 




Sp. Type 


Ref. 


Multiplicity 


WR 


11 


W(J8 + 08.5 


a 


Confirmed binary with a period of 78.5 d 


WR 


14 


WC7 


b 


Suspected binary 


WR 


21a 


WN6 + 03 


c 


Confirmed binary 


WR 


39 


WC7 


b 


Possible visual binary 


WR 


48 


WC6 + 09.5 


a 


Confirmed binary; suspected triple system 


WR 


79a 


WN9ha + ? 


d 


Visual binary 


WR 


89 


WN8ha + OB 


d 


Visual binary 


WR 


90 


WC7 


b 


No evidence for binarity 


WR 


98a 


WR + OB 


c 


Confirmed binary with a period of 










about 1 yr 


WR 


104 


WR + OB 


c 


Confirmed binary with a period of 










about 1 yr 


WR 


105 


WN9 


a 


No evidence for binarity 


WR 


112 


WC9 


a 


No evidence for binarity 


WR 


125 


WC7 + 09 


a 


Confirmed binary with a period of > 15 yr 


WR 


137 


WC7 + OB 


a 


Confirmed binary with a period of > 13 yr 


WR 


140 


WC7 + 05 


a 


Confirmed binary with a period of 










about 7.9 yr 


WR 


146 


WC5 + 08 


a 


Visual binary 


WR 


147 


WN8 + B0.5 


a 


Visual binary 
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In the case of B-type stars, only a few objects have been proposed to be 
non-thermal radio emitters. However, the non-thermal nature of at least a 
fraction of their radio emission is uncertain and still needs to be confirmed. 
Considering, the lack of relevant observational indications for a non-thermal 
radio emission from B stars, this category will not be considered here. 

The issue of the binarity of many of the objects listed in Tables[T]and[2]has 
been elucidated thanks to spectroscopic observations in the optical domain. 
In some cases, the multiplicity was studied following high angular resolution 
techniques. However, our knowledge of the multiplicity of these stars is not 
complete. In some cases, the multiplicity has not yet been investigated, and 
in several cases the presence of the companion has only been detected and a 
lot of information is still lacking concerning the members of these multiple 



Table 2 List of the non-thermal radio emitting O-type stars. The information 
on the non-thermal nature of the radio emiss i on fro m th e targets collected i n 
this table are taken from (a) [Contrcras et al (1997]), ( b) IBieging et al (11 989), 
(cl ISetia Gunawan et"al (I2003T). (d) IDrakei (11990). ( el IBenaglia and Korib alski 



2004), (f) Leitherer et all ll 19951). ( g) IBenaglia et all (2001a), (h) IBlomme et all 
(2007) and (i) IBlomme et all (|2005f ). Fo r the multiplicity and/ or the spectral 
types, the i nformation are take n from (\) IRauw et a l (1999), fk) IDe Becker et all 
(120041 JS. (T IRauw et all H2005blh (m) IGies et all (119931). (n) IHarvin et all (120021) . 
(o)IDe Becker et all (12006al ) . (p)INelan et al 112001 ). (q) IRauw et all (120031). (r)IGiesl 
(119871 ). (s) Gossct (2005) (t) Leith erer etaM l987). and (u) IDe Becker et all (1200 l;il). 



Star 


Sp. Type 


Hef. 


Multiplicity 


CygOB2#5 


06-7Ib + Ofpe/WN9 


a 


Confirmed binary, probably 
triple 3 ' • a 


CygOB2#8A 


06If + 05.5111(f) 


b 


Confirmed binary with a 
period of 21.908 d fc 


CygOB2#9 


05If 


b 


Multiplicity currently 
under investigation 


CygOB2-335 


07V 


c 


Multiplicity not yet 




04f+ + ? 




investigated 


9Sgr 


b 


Confirmed binary with a 
period of a few years' 


15Mon 


07Ve + ? 


d 


Confirmed binary" 1 


50riA 


09.511 + ? 


b 


Confirmed multiple system™ 


aOriAB 


09.5V + ? 


d 


Confirmed multiple system d 


HD 15558 


05III + 07V 


b 


Confirmed binary with 
a period of about 440 d, 
possibly triple 


HD 93129A 


02If + ? 


c 


Visual binary p 


HD 93250 


03V((f)) 


f 


Suspected binary" 


HD 124314 


06V(n)((f)) + ? 


g 


Confirmed binary r 


HD 150136 


05IIIn(f) + ? 


g 


Confirmed binary, probably 
triple" 


HD 167971 


05-8V + 05-8V + 081 


b,h 


Confirmed binary, 
probably triple 








HD 168112 


05III(f+) 


b,i 


Suspected binary" 


CD-47°4551 


05If 


g 


Multiplicity not yet 
investigated 
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systems. One of the main problems relevant to the issue of the multiplicity of 
stars in general is the strong bias that affects the census of known binaries. 
This problem is twofold: (1) A first origin of this observational bias is the 
problem of the orbital period: typical observing runs last a few days, and are 
often separated by several months or years. A consequence is that most of the 
known binaries have periods of a few days, and a few others are just referred 
to as long period binary candidates because long term radial velocity shifts 
were detected on time-scales of one or several years. That is the reason why 
we have only little information on the multiplicity of several targets discussed 
here. It appears that intensive spectroscopic monitoring should be performed 
to investigate poorly explored regions of the orbital parameter space, in or- 
der to probe some of the fundamental properties of massive stars. Moreover, 
some of the systems quoted in Table[2] appear to be multiple systems whose 
spectroscopic investigation is not straightforward, pointing out the necessity 
to develop or improve spectral disentangling tools. (2) A second origin of 
this observational bias is the inclination of the systems we are talking about. 
As most of the binaries are revealed by spectroscopic monitorings aiming at 
the measurement of radial velocity excursions, systems characterized by low 
inclinations are often erroneously classified as single stars. In this context, a 
promising technique likely to provide a wealth of information is interferom- 
etry. Thanks to fairly long baseline interferometers such as the Very Large 
Telescope Interferometer (VLTI, in Cerro Paranal, Chile) with near infrared 
spectro-imaging capability, we may for instance confirm the binary status 
of some stars, and study the orbital motion of long period systems whose 
orbital parameters have not yet been investigated. Considering the current 
angular resolution achievable with such interferometers (of the order of a few 
mas), intcrfcromctry should be devoted to multiple systems characterized by 
rather large separations (at a distance of 1 kpc, 1 AU of linear separation 
corresponds to only about 1 mas) . 

The information that can be obtained from such multiplicity studies con- 
cern (1) the geometry of the binary or higher multiplicity systems and (2) 
the spectral and wind properties of their members. Both are crucial in order 
to constrain the properties of the non-thermal emission components likely to 
be studied in the radio and high-energy domains (see Sect. [6]). 



Finally, I should emphasize that the binary/multiple scenario is a valu- 
able working hypothesis which agrees with almost all of the observational 
constraints available at the time of the writing of this review. More ob- 
servations are needed to validate - or invalidate - this hypothesis. These 
additional data are mostly required in the case of O-type stars. Indeed, two 
members of the O-type non-thermal radio emitters catalogue presented in Ta- 
ble[5]are only suspected binaries (i.e. HD 168112 and HD 93250), and recent 
spectroscopi c observations failed to r eveal the presen ce of a companion (see 
respectively iRauw et all (|2005bl ) and iLanottd ( 20061 )). In other cases, such 
a spectroscopic monitoring is still lacking. If these stars - and other non- 
thermal radio emitters - are indeed single stars, an alternative scenario will 
be required to complement the one that is described here. However, as the 
general scenario described in the previous sections in the context of binary 
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systems still seems to be the most appropriate, such an alternative scenario 
will not be considered here. 



8 Multiwavelength study of the non-thermal emission 

Considering the detailed discussion above, we see that a lot of physical aspects 
should be taken into account simultaneously in order to address the issue of 
the non-thermal emission from massive binaries in various domains of the 
electromagnetic spectrum. For this reason, a multiwavelength approach is 
definitely needed. In this section, I describe several results obtained in various 
energy domains, and I discuss some prospects for future observations. 



8.1 The radio domain 



As mentioned in the introduction, the non-thermal emission from massive 
stars was first reported in the radio domain. Many radio observations con- 
tributed to the census of non-thermal radio emitters that was made by the 
end o f the 80's (IWright and Barlowll975l : lAbbott et a]ll984ll986t lB ieging et all 



I1989D . In the large radio survey published bv lBieging et aT(jl989l ). about 25 % 
of the early-type stars detected in the radio domain presented the signature 
of a non-thermal radio spectrum. Since then, this fraction has increased up 
to a value close to 40 %. The present census of non-thermal radio emitters is 
described in Tables [1] and [2] respectively for WR and O-type stars. 

Beside multifrequency radio observations of massive stars able to re- 
veal flux densities and spectral indices, crucial information came from high- 
angular resolution radio observations (VLA, VLB A, MERLIN). Thanks to 
these data, major progres has been made by resolving non-thermal emis- 
sion regions in the case of a few colliding-wind binaries, therefore confirming 
that the non-thermal emission arises from the wind-wind interaction region 
and not from the stars themselves. Among th ese examples, we may consider 
CvgOB2#5 (IContreras et allll997h WR146 (lO'Connor et alll20ol . WR147 
([Williams et a]|ll997f ) andWR 140 (|Doughertv et alll2005ah . 



An important point worth mentioning here concerns the confrontation 
of observational rad io results with detai led simulations of radio emission 
from single O-stars. Van Loo et a 3 (|2006l) demonstrated that the observed 



non-thermal radio emission from O-stars cannot be explained by current 
models of wind-embedded shocks produced by the instability of the line- 
driving mechanism, therefore lending strong support to the binary scenario 
for the non-thermal radio emission from massive stars. This latter scenario is 
thus supported by both theoretical (radio emission models) and observational 
(multiplicity studies) arguments. 
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8.2 The soft X-ray domain 

One of the main questions that has been addressed in the context of ongo- 
ing multiwavelength observations is that of the high-energy counterpart to 
the non-thermal emission already detected in the radio domain. As shown 
in Sect. El several processes may be responsible for the production of a non- 
thermal high-energy emission. Among the processes discussed in Sect. El in- 
verse Compton scattering is expected to be the most significant one as it 
requires a large amount of UV and visible photons, and massive stars are 
strong sources over precisely these energy domains. As a consequence, we 
may expect a non-thermal component in the high-energy domain following a 
power law spectral shape: 

F(E) oc E a (35) 

The spectral index (a) is directly related to the index of the relativistic 
electrons: 

a = ^ (36) 

As a result, in the case of strong shocks where n is equal to 2 (see Eq.EJ, 
we expect the non-thermal emission component in the high-energy domain 
to have a spectral index of 0.5. 

Several non-thermal radio emitting massive stars have been observed with 
the European X-ray satellite XMM-Newton that is sensitive at energies be- 
tween 0.5 and 10.0 keV. As mentioned in Sect.JU massive stars and massive 
binaries are thermal emitters in the X-ray domain. These thermal emission 
components are characterized by plasma temperature s intimately depen dent 
on the pre-shock velocities of the plasmas that collide ( Stevens et allll992l) . In 



the case of intrinsic shocks in individual stellar winds, the pre-shock velocities 
are of the order of a few 100 km s -1 . The post-shock temperature is therefore 
of the order of a few 10 6 K. In the case of colliding winds, the situation is a 
bit more complicated. If we consider long period binaries, i.e. binaries with 
orbital periods longer than a few weeks, the separation between the stars in 
the binary system is rather large. As a consequence, the winds have enough 
space (and time) to accelerate up to their asymptotic velocity (the terminal 
velocity) before colliding. In such a situation, the typical pre-shock velocity 
is of the order of 2000 km s -1 and the post-shock temperature is of the or- 
der of several 10 7 K. In the case of close binaries, i.e. binaries with periods 
not longer than a few days, the winds collide in their acceleration zone. The 
pre-shock velocities are therefore much lower than in the case of long period 
binaries, i.e. of the order of 1000 km s -1 , and the post-shock temperature is 
of the order of 10 7 K. 

As thermal processes dominate the X-ray spectrum of massive stars in the 
XMM-Newton bandpass, we may expect the thermal emission components 
to have a strong impact on our capability to detect a non-thermal emission 
component in the X-ray spectrum of massive binaries: a putative non-thermal 
emission component will be detected only if it is not overwhelmed by the ther- 
mal emission. 
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Several O-t ype non-thermal radio emitter s have been observed wi th XMM- 
Newton: 9 Sgr feauw et al2002h . HD 16 8112 (iDe Becker et all2004al) . HD 167971 
(|De Becker et all2005al) . CvgQB2 #8A (IDe Becker et all2006bl) and Cyg OB2 #9 
(jRauw et al 2005a). These targets are confirmed or suspected long period bi- 
naries. In each case, we observe an X-ray spectrum characteristic of rather 
high plasma temperatures. Moreover, the spectrum is hard, i.e. the X-ray 
flux descreases rather smoothly as a function of the energy. None of these 
targets revealed unambigously a non-thermal X-ray contribution in their X- 
ray spectru m below 10.0 keV. I n addition, it is worth mentioning that the 
analysis by I Skinner et all ([2007) of XMM-Newton data of the very long pe- 
riod WR + O binary WR147 reveals a thermal spectrum below lOkeV. On 
the other hand, the XMM- Newton observation of the short period colliding- 
wind binary HD 159176 (|De Becker et alll2004ct ) revealed a rather soft ther- 
mal spectrum along with a power law suggesting the presence of a non- 
thermal emission component. If it is confirmed, HD 159176 would be the 
first O-type non-thermal emitter in the X-ray domain below 10.0 keV. This 
result is therefore very interesting, mostly if we consider that this latter tar- 
get is not a non-thermal radio emitter. A second candidate for non-thermal 
emission from a short period binar y in the soft X-ray domain is FO 15 
([Albacete Colombo and Micelall2005T) . 



Considering these observational results, we may distinguish two main 
categories. The first category is that of long period binaries, i.e. binaries with 
periods longer than a few weeks. In this case, the separation is large and the 
stellar winds reach their terminal velocities before they collide. The post- 
shock temperature is high and the thermal emission is hard. In addition, 
the collision zone is far from the stars and the UV radiation field is weak 
in the interaction zone. This situation is therefore not favourable to the in- 
verse Compton scattering process that requires a large amount of low energy 
photons to be efficient. The putative non-thermal emission is thus expected 
to be weak. As a result, long period binaries are unlikely to reveal any non- 
thermal contribution in their soft X-ray spectrum (see Fig. 02 part (b)). In 
the radio domain, the large binary separation allows a significant fraction of 
the synchrotron radio emission to emerge from the combined stellar winds, 
and such systems may appear as non-thermal radio emitters. 

In the case of short period binaries, i.e. binaries with periods not longer 
than a few days, the binary separation is shorter. As a consequence, the 
post-shock temperature at the collision zone is lower and the thermal X-ray 
spectrum is softer, i.e. the X-ray flux decreases more rapidly as the energy in- 
creases. Moreover, the short distance between the stars and the collision zone 
is in favour of the inverse Compton scattering mechanism. The relative in- 
tensity of the putative non-thermal emission component may be rather high. 
In this situation, a non-thermal emission component may be revealed in the 
X-ray domain (see Fig. [5) part (a)). In the radio domain, the synchrotron 
emission is unlikely to be unveiled as the emission region is strongly burried 
in the combined stellar winds, and therefore significantly absorbed. As a re- 
sult, short period binaries will not be detected as non-thermal radio emitters. 
Such systems may be referred to as synchrotron- quiet non-thermal emitters, 
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Fig. 5 Effect of the hardness of the thermal X-ray spectrum on the detectability of 
a putative non-thermal emission component. Cases (a) and (b) stand respectively 
for short period and long period binaries. On the left: schematic view of the binary 
system, with the bold double arrows emphasizing the significant difference in the 
separation between the star and the collision zone (CZ) depending on the case (long 
or short period). On the right: simulated XMM-Newton spectrum between 0.5 and 
10.0 keV including thermal (Ti and T2) and non-thermal (NT) components. The 
temperature of the Ti component and the spectral index of the NT component 
are the same in the two cases. The only difference resides in the temperature of 
the second thermal component, which is higher in case (b) than in case (a). In the 
simulated spectrum of case (b), the strong spectral line at about 6.7 keV is due 
to Fe, and is often considered as a signature of a high plasma temperature. Case 
(a): the power law is revealed, at least above about 3keV. Case (b): the power law 
is overwhelmed by the high temperature thermal component, and is therefore not 
revealed. 



provided they produce a detectable non-thermal radiation in the high energy 
domain but not at radio wavelengths. 

In summary, it appears that the simultaneous detection of a non-thermal 
emission in the radio domain and in the X-ray domain below 10.0 keV is 
unlikely. However, as the detection in the X-ray domain is mostly inhib- 
ited by the presence of the strong thermal emission, a simultaneous detec- 
tion can be envisaged if we investigate the higher energies - the hard X-rays 
(above 10.0 keV) and the soft "/-rays - where the thermal emission is absent 
(jDe Beckerll2005h . 
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8.3 The hard X-ray and soft 7-ray domains 

The detection of a high-energy counterpart to the non-thermal radio emis- 
sion still needs to be achieved in the ca se of O-type stars. The analysis 
of INTEGRAL-lSGRl data performed by iDe Becker et~a1 (|2007bl) between 



20 keV and 1 MeV failed to reveal any hard X-ray of soft 7-ray emission 
related to the O-type star population of the Cyg OB2 region. 

Another point worth mentioning concerns the magnetic field. As shown 
in this paper, the magnetic field intervenes in differents aspects of the ques- 
tion of the non-thermal emission from massive binaries (synchrotron radio 
emission, acceleration of particles...). One can therefore wonder whether the 
study of the non-thermal emission from massive stars could lead to an indi- 
rect method to estimate the magnetic field strength, at least locally, in the 
collision zone. Indeed, when we consider for instance Eq. ll4l we see that in 
the case of a binary system, the simultaneous knowledge of several parame- 
ters such as the bolometric luminosities, the distance between the stars and 
the interaction zone, the synchrotron radio luminosity and the high-energy 
inverse Compton luminosity, may allow to estimate the strength of the local 
magnetic field. However, as discussed above, the simultaneous determination 
of the latter two quantities requires the combination of radio and high-energy 
(hard X-rays and soft 7-rays) observations. This prospect underlines the in- 
terest of observing non-thermal radio emitters at higher energies. 



8.4 The very high-energy 7-rays 

When we are dealing with very high energy processes, we should keep in 
mind the fact that energy of the particles involved in these processes must 
be very high as well (i.e. particles can only transfer to radiation the energy 
they possess). As discussed in Sect. 16 .51 the maximum Lorentz factor that can 
be reached by relat ivistic electrons in colliding- wind bi naries should be of the 
order of 10 4 to 10 5 (jDe Beckerl200l iPittard et alll2006f) . As a result, the high- 



energy photons produced by electrons through inverse Compton scattering 
should not reach energies higher than a few GeV. However, we should not 
forget that protons are also accelerated. In Eq.[32j we see that the maximum 
Lorentz factor reached by particles depends on their mean free path in the 
acceleration zone, i.e. the typical length scale between two scattering centers. 
Generally, it is considered that this mean free path is directly related to 
the gyroradius of the particle. The more massive the particles (protons are 
about 1800 times more massive than electrons), the shorter the gyroradius 
and therefore shorter is their mean free path. The confinement of protons in 
the acceleration zone will therefore be significantly improved as compared to 
that of electrons. Protons are thus expected to reach much higher energies 
than electrons. 

Considering th e depend e nce o f the Thomson cross section on the mass 
of the particle (see iLongairl ( 1992T )). inverse Compton scattering by protons 
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will be inefficient. As a result, other hadronic processes^ such as neutral pion 
decay should be considered (see Sect. [5]). This latter proces s has indeed been 
proposed to be at work in the vicinity of mas sive stars ( Chen et"al Il996t 
iReimer et all2006t iPittard and Doughertvll2006l ). 

With the advent of ground-based Cherenkov telescopes able to detect 7- 
ray photons with energies between about 100 GeV and a few TeV, a new 
window of high-energy astrophysics is now open. Among the so-called TeV 
sources that have been detected in the recent years, some of them may be 
related to early- type star populations. The first unidentified TeV source was 
indeed discovered in the rich OB as sociation CygOB2 by the High Energy 
Gamma Ray Astronomy consortium ( Aharonian et a3l2005l ). and its position 
has been slightly revised re cently thanks to obser vations with the Whipple 
Observatory 10 m telescope ([Konopelko et alll2007l) . The nature of this source 
is still unknown, even though it has been proposed that it may be related to 
the rich populatio n of massive stars in the CygOB2 region ( Butt et al 2003: 
iTorres et all 12004 ). It is worth noting that the recent study by iButt et all 
(|2007 ) revealed the existence of a radio counterpart to this TeV source. 
Another TeV so urce was discovered cl ose to the Galactic Center using the 
same technique (I Aharonian et alll2004j ). This latter source was proposed by 
lOuataert and Loebl ((20051 ) to be related to the massive star population lo- 
cated in the central region of the Galaxy. However, a detailed study of the 
massive star population close to the Galactic Center is still lacking, although 
some recent investigations show ed that tens of early- type stars (O-type and 
WR) are present in that region (jPaumard et a 3 12004 ). 



9 Recent theoretical developments 

The many observational results gathered concerning the non-thermal radio 
emission from massive binaries, along with the prospects for the detection 
of a high-energy counterpart, motivated astrophysicists to develop models 
devoted to this issue. Significant advances have been made during the last 

years. 

In the radio domain, recent models bv lDoughertv et all (|2003l l2005bh im- 



proved significantly our vision of the physical processes at work in colliding- 
wind binaries. These new approaches use hydrodynamical codes to describe 
the stellar winds and the wind collision regions in order to model the syn- 
chrotron emission, taking into account several crucial effects such as the Razin 
effect or the fre e-free absorpt i on by the wind material. More recently, the de- 
tailed work by IPittard et all ( 2006t ) improved significantly the modelling of 



non-thermal emitting colliding winds through the inclusion of the inverse 
Compton cooling. The strong impact of this latter process is indeed likely to 
alter substantially the properties of the population of relativistic electrons, 
and therefore the synchrotron emission as well. These models predict also a 
steepening of the electron spectrum, i.e. an electron index larger than 2 (see 

5 Hadronic processes are processes involving hadrons (i.e. protons, neutrons or 
ions) , in contrast to leptonic processes requiring the intervention of leptons such as 
electrons. 
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Sect. [3]) , where IC cooling is strong. This comes from the fact that IC cooling 
appears to be more efficient at cooling the higher energy electrons. It is inter- 
esting to note that this may explain why the spectral index detected in the 
case of HD 159176 - the first O-type binary possibly presenting a non-thermal 
component in the soft X-rays - is significantly larger than 0.5 which is the ex- 
pected theoretical value for strong shocks. The inclusion of IC cooling allows 
moreover to estimate the non-thermal high-energy spectrum of colliding- wind 
binaries following a more satisfactory approach than that allowed by Eq.[T4l 
These recent studies allow indeed to fit models to radio data in order to deter- 
mine the population and spatial distribution of relativistic electrons without 
an a priori knowledge of the local magnetic field. In addition, it appears that 
the relativistic particles themselves are expected to alter the hydrodynamic 
properties of the shocks responsible f or the particle acceleratio n, like in the 
case of supernova remnants (see e.g. iDecourchelle et a3 (2000)). When this 
shock modification is efficient, the relativistic electron population is characte r- 
ized by an index (se e Eq.[H) larger that 2 (se e lPittard and Dougherty! ( 20061 )). 
On the other hand. IPittard and Dougherty] (|2006l ) suggested that in the case 
of WR 140 (see below), the energy spectrum of the relativistic electrons may 
be affected by the multiple accelerations likely to occur consecutively in dif- 
ferent shocks. In this scenario, the electrons are first accelerated by intrinsic 
shocks in individual stellar winds, before being further accelerated by the 
shocks of the wind-wind interaction region. When this re- acceleration pro- 
cess occurs, the spectral inde x of the relativistic electrons is lower than 2 
( Pittard and Doughertvl[2006l ). 

Alternative models have also been developed by iReimer et all ( 20061 ) in 
order to address the issue of the non-thermal high-energy emission from 
colliding-wind binaries. This latter approach includes the particle accelera- 
tion self-consistently and allows to compute high-energy spectra through var- 
ious emission processes such as inverse Compton scattering (in the Thomson 
and Klein-Nishina regimes) or neutral pion decay. Such models are necessary 
in order to prepare future high-energy observations with space-borne obser- 
vatories such as GLAST, or ground-based Imaging Atmospheric Cherenkov 
Telescopes (IACTs). 



10 Two particularly interesting systems: WR 140 and 
CygOB2#8A 

In this section, I discuss more specifically the case of two non-thermal radio 
emitters that deserve particular attention. These two cases have been selected 
to span a rather large domain of the stellar and orbital parameter space. The 
first one is a long period binary system made of a WR and of an O-type star 
(WR140), and the second one is a shorter period binary constituted of two 
O-type stars (CygOB2#8A). 
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10.1 WR140 



This interesti ng binary was the first star for which non-thermal emission 
was detected (jFlorkowski and Gottesmanlll977h . and it has long been con- 
sidered as the archetypal non-thermal emitting WR star. It consists in a WC7 
primary with an 04-5 secondar y orbiting with a per iod of about 7.9 yr in a 
highly eccentric orbit (e = 0. 88). I Williams et all (jl990h showed that the highly 



elliptical orbit is responsible for strong variations all the way from X-rays to 
radio wav elengths. The intens i ve lon g term radio monitoring of WR 140 per- 
formed bv lWhite and Beckerl (199.5) provided for the first time a radio light 



curve of a non-thermal radio emitter during a complete orbital cycle. This 
light curve presented strong modulations unlikely to be attributed to free- free 
absorption only. The variations are believed to be at least partly due to in- 
trinsic changes in the physical conditions underlying the intri nsic synchrotron 
emissio n as a function of the orbital phase. More recently, iDoughertv et all 
(|2005ah provided a more detailed view of the dependence of the properties of 



the synchrotron radio spectrum of WR 140 as a function of the orbital phase 
through multifrcquency observations. In addition, the latter authors reported 
on an arc of emission - reminiscent of the bow-shaped morphology expected 
for the wind collision region - that is rotating as the orbit progresses, even 
though it should be noted that the observed morphology is also affected by 
the absorption of radiation by the wind material close to the shock. 

WR 140 has been proposed as a possi ble counterpart to the unide ntified 
EGRET 7-ray source 3EG J2022+4317 (|Benaglia and Romero! f2003h . even 



though the binary lies at the limit of its error box. The investigation of a large 
set of INTEGRAL-IBIS data did not lead to the d etection of a hard X-ra y 
emission (between 20keV and 1 MeV) from WR 140 (|Pe Becker et alll2007bD 
On the basi s of models such as discussed in Sect.l9| IPittard and Doughertvl 
( 20051 l2006h calculated the expected high-energy spectrum of WR 140. From 



the fit of the radio data at a particular orbital phase, these authors estimated 
that the high-energy emission would be dominated by IC scattering up to 
50 GcV, while the emission from hadronic processes may reach a few TeV. The 
predicted emission levels suggest promising results from future observations 
with GLAST and Cherenkov air shower telescope arrays. A future detection 
of WR 14 with these latter instruments is also strongly suggested by the 
results of iReimer et a 3 l|200fift . 



10.2 CygOB2#8A 



This particularly interesti ng target was unam biguously classified as a non- 
thermal radio emitter b y iBieging et all ( 19891 ). The spectroscopic study of 
iDe Becker et "all fl2004ri ' led to the discovery of its binar i ty, w ith an or- 
bital period of about 22 days, and IDe Becker and Rauwl ( 2007f ) provided 
some indications of a wind-wind interaction in this system. In the X -rays 
below lOkeV, the stud ies performed for instance by IDe Beckerl ( 200l[ ) and 
IDe Becker et all ( 2006bf ) did not reveal any non-thermal emission component. 
The soft X-ray spectrum appears indeed to be essentially dominated by the 
strong thermal and phase-locked variable emission from the colliding-winds 
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( De Becker and Rauwl 1200,4 iDe Beck er et a jl l2006bh . In the radio domain, 
the recent monitoring performed bv lBlommd ( 2005h revealed a strong phase- 
locked variab ility of the radio flux, in excellent agreement with the ephemeris 
published bv lDe Becker et al (|2004bl h As the radius of the sphere with unit 
optical depth in the radio domain is expected to be very large (of the or- 
der of 100R*jl, it was a priori unexpected to detect the synchrotron radio 
emission that is most probably produced in the wind interaction region, i.e. 
deep inside the combined winds. The fact that we clearly detect it is most 
probably due to a combination of wind porosity and inclination effects. A 
more detailed monitoring of the orbital cycle in the radio domain is planned 
in the future in order to constrain the physical cirucumstances underlying 
the non-thermal emission processes at work in the colliding-wind region. 

A point that deserves to be mentioned is that CygOB2^=8A - along with 
other non-thermal radio emitting massive stars - lies inside the error box of 
the unidentified EGRET source 3EG 2032+4118. This 7-ray source has al- 
ready been proposed to be associated to m assive stars in th e CygOB2 asso- 
ciation: CygOB2#5 (jBenaglia et alll2001bl ). CygOB2#8A (|De Becker et all 
2005b). However, recent observations with the IBIS imager onboard the Euro- 
pean INTEGRAL sa tellite between 20 keV a nd 1 MeV did not allow to detect 
the EGRET source (De Becker et a 1 l2007bh . Observations with future hard 
X-ray observatories benefitting from a better sensitivity than INTEGRAL 
are required to address the issue of the non-thermal high-energy emission 
from this promising system. In this context, observations between 10 and 
100 ke V with the SIMBOL-X satellite are very promising ( De Becker et all 
l2007ah . 

Among the O-type non-thermal radio emitters, Cyg OB2 #8 A appears to 
be the best candidate for an application of state-of-the-art models (jDoughertv et all 
120031; iPittard et alll2006l ). Cyg OB2 #8 A is indeed the O-type non-thermal ra- 
dio emitter with the best constrained st ellar, wind, and orbital parameters 
( De Becker et a 1 l2004bl: IDe Beckerl l20"05h . and constitutes therefore a priv- 
ileged target to probe the physics underlying the non-thermal processes in 
massive colliding-wind binaries. 



11 Concluding remarks 

Since the discovery of the occurrence of non-thermal processes in several 
massive stars, significant progresses have been achieved. First, the nature of 
non-thermal radio emission seems to be established: it has been identified 
as synchrotron radiation. Second, we are more and more convinced that the 
binary scenario is the most adapted to explain the non-thermal emission 
observed in the radio domain. 

A crucial point that is extensively demonstrated in this review is the 
necessity to address the issue of the non-thermal emission from massive bi- 
naries following a multiwavclcngth approach. The discussion of the two best 

6 The radius of the sphere with unit radio optical depth is a decreasing function 
of the frequency. It may be of the order of a few tens of stellar radii at 3 cm, and 
it can reach several hundreds of stellar radii at 20 cm. 
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studied members of the catalogues presented in Tables Q] and [2] - WR 140 
and Cyg OB2 #8A - shows that many types of information on these systems 
(stellar parameters, wind properties, orbital elements) are required in order 
to discuss in detail their non-thermal emission in several energy domains, 
through the application of recent theoretical models. These models are now 
being applied to a few WR-type non-thermal radio emitters. The next step 
will consist in their applicaiton to O-type systems such as CygOB2#8A. 
Such an extension to O-type binaries is a necessary condition to probe a 
more extended domain of the parameter space. In addition, the investiga- 
tion of synchrotron-quiet massive binaries should not be neglected as the 
latter may also constitute non-thermal radiation sources in the high-energy 
domain. As a result, it may be more justified to talk about a more extended 
class of objects called non-thermal emitters, among which one could find sys- 
tems detected as non-thermal radio sources - the so-called non-thermal radio 
emitters - and possibly non-thermal emitters in the high-energy domain that 
are not detected as synchrotron emitters in the radio domain. 

Even though leptonic processes are by far the most studied in this con- 
text, the prospect of very high 7-ray emission from massive stars deserves 
to consider hadronic processes as well. State-of-the-art models now include 
both processes simultaneously, and provide powerful tools likely to sketch a 
more complete view of the multiwavelength spectrum of non-thermal emit- 
ting early-type colliding-wind binaries. Additional studies aiming at the im- 
provement of the knowledge of useful parameters are required, notably in 
the visible and near-infrared domains. Such studies are needed in order to 
investigate for instance the acceleration of particles in colliding winds of early- 
type stars, and whether massive binaries may contribute to very high-energy 
sources such as the TeV emitters that have been detected u sing Cherenkov 
telescopes in the Cyg OB2 association (|Aharonian et aj|2005h or close to the 
Galactic Center ([Aharonian et a]|l2004l ). 
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